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Abstract: For very sparse tree land individual tree was the basic element of interaction between atmosphere and 
the surface. Drag of isolated tree was preliminary aerodynamic index for analyzing the atmospheric boundary 
layer of this kind of surface. A simple pendulum method was designed and carried out in wind tunnel to measure 
drag of isolated tree models according to balance law of moment of force. The method was easy to conduct and 
with small error. The results showed that the drag and drag coefficient of isolated tree increased with decreasing of 
its permeability or porosity. Relationship between drag coefficient and permeability of isolated tree empirically was 
expressed by quadric curve. 
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Introduction 

Research on the interaction between vegetative 
surface and atmosphere is one of the aims of at¬ 
mospheric boundary layer meteorology. Momentum 
transfer between the atmosphere and vegetative 
structure is a function of the physical attributes of the 
structure and the wind field (Grant 1985) drag of 
vegetation is an important index of the function. 

There are various vegetative structures in surface, 
natural or planted or mix of them. Homogeneously 
dense vegetation may be the simplest surface in 
structure. Its drag was studied in many papers 
(e.g.Thom 1971; Raupach 1979; Bosveld 1997). An¬ 
other structure of frequently researched was 
shelterbelt or windbreak, as it was widely constructed 
to reduce soil evaporation or wind speed in arid or 
windy areas, Its drag also was evaluated in many 
studies (e.g. Hagen 1971; Seginer 1975; Takahashi 
1978; Wilson 1997; Zhou 1994; Zhu 1992; Guan 
1995, 1996, 1998). Homogeneously dense vegeta¬ 
tion was considered entirely as rough surface and 
shelterbelt or windbreak was looked as porous ob¬ 
stacle belt. Feature of individual tree was not taken 
into account in the studies. 

Scattered or sparse trees is often seen a kind of 
surface, such as savannah, agroforest, parks and so 
on. Individual trees can be taken as elements of in- 
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teraction between atmosphere and surface. Evalu¬ 
ated theoretical models of boundary layer parameters 
of this kind of surface (Raupach 1992, 1995) need 
the drag of isolated tree as an input index. Some 
previous studies dealt with individual tree (e.g. May- 
head 1973) or organs (e.g. Grant 1985; Vogel 1989; 
Holland 1991). But the measurements in the studies 
had few replications, thus permitted no reasonable 
extrapolation of the results to other kind trees. 

This study is an attempt at measuring drag of iso¬ 
lated tree model in wind tunnel and describes the 
relationship between drag of isolated tree and its 
structural index. 

Optical porosity and aerodynamic porosity of 
isolated tree 

Two parameters, optical porosity and aerodynamic 
porosity, are introduced to describe the geometric 
characteristic and aerodynamics of isolated tree. The 
definition of the parameters and relationship between 
them were expounded in another paper (Guan and 
Zhu 2000). 

Wind tunnel experiments on drag of isolated 
tree models 

Wind tunnel and tree models 

A simple wind tunnel of blow down type was made 
for the experiments. The length of test section is 10.0 
m, width is 0.8 m, and height is 0.6 m. A honeycomb 
and 4 wire screens were installed at the upstream 
tunnel entrance to diminish flow rotation and large 
scale eddying. It is capable of generating a maximum 
air velocity of 8.0 m/s in the test section. Wind veloc- 
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ity probe (probe of hot bulb anemometer) inserted 
through the holes open in tunnel ceiling and changed 
the height of measuring position. Three spires (40 cm 
height) were placed at upstream and grits (average 
diameter 3 mm) were uniformly spread on tunnel floor 
to product thick boundary layer. Measurements were 
carried out at the test section near exit where the 
boundary layer thickness was 14-16 cm and the 
roughness length was 0.018 mm (S.D.=0.006 mm). 

A tree model was consisted of crown elements 
(plastic bunches) and trunk (iron wire), which string 
the element. Three kinds of elements were used to 
construct tree models. The crown became dense or 
porous as the elements increase or decrease. All 
models were designed in 8 cm height. 10 kinds of 
trees were combined with different figures. Solid tree 
model was substituted by wood cylinder. Structural 
characteristics were described in another paper 
(Guan and Zhu 2000). 

The experiment description 

Because no instruments were available for meas¬ 
uring so small drag of the tree models, the pendulum 
method was designed to the experiment (Fig.1). The 
ceiling at the exit of the wind tunnel was opened a 
long hole (1 cmx 10 cm) along the wind direction. A 
pendulum was designed whose rigid arm was freely 
suspended just into the hole and didn’t touch the 
ceiling board. The parallel double hooks at upper end 
of the arm hooked on a triangle section beam hori¬ 
zontally were fixed at steel frame above tunnel ceiling. 
And the up ridge of the fixed beam was sharply pol¬ 
ished to decrease the friction with the hook as the 
pendulum vacillated. The double hooks controlled the 
arm only moving along the hole (wind direction). The 
down end of the arm is apart from the wind tunnel 
floor about 8 cm (height of the tree models) when it 
freely suspended. On tunnel floor just down the arm a 
ruler was drawn. A steel needle fixed at the down end 
of the arm as indicator. When the arm vacillated with 
the wind the indicator indicated the displacement of 
the arm end. The arm moved back and forth even in 
uniform wind. So a detention device installed out of 
the wind tunnel to detain the vibration of the arm. The 
device extended a detention arm into the tunnel to 
stop the pendulum arm upward. 

In the experiment the position of the end of pen¬ 
dulum arm was read first in still air. Then the uni- 
form-velocity wind flow was generated and detention 
arm was gradually regulated out of the tunnel to the 
optimum position that the pendulum arm was slightly 
bumping on the detention arm. The displacement of 
the end of pendulum arm was read. After stopping 
the flow the position of the end of pendulum arm in 
still air was read again to test the pendulum arm 
whether back to the initial position. 


Then the needle was detached from the arm and 
fixed to the trunk base of the tree model also as in¬ 
dicator of displacement. And the tree model fixed to 
down end of the pendulum arm and maintained the 
trunk base departing from the tunnel floor about 1mm. 
The measuring process repeatedly described above 
and displacement was recorded. The atmospheric 
pressure and air temperature also were recorded as 
the experiment going on. The weight of pendulum 
arm, tree models and the needle were weighing. 

In order to control the displacement of pendulum 
arm end in small scale, the more weight pendulum 
arm was selected, as the wind was stranger. This 
could diminish the errors that generated from much 
dip and elevation of the tree model. In the experi¬ 
ments the length of pendulum arm was 160 cm and 
the displacement of pendulum arm end was less than 
7.2 cm, so the swing angle was less than 2.6° 

Calculation of the drag of tree model 

The force diagram of the pendulum system was 
also shown in Fig.1. When no tree model attached, 
the pendulum arm was forced by 

(1) Gravity G L . Its center was at the middle point O 
of the arm. 

(2) Drag F L . Its center was at the point M L of the 
arm. And the dip angle was cx L and down end dis¬ 
placement was di_. 

The diameter of needle was 0.2 mm, only 1/40 of 
that of the pendulum arm. So the gravity and drag of 
the needle were neglected as they much smaller than 
that of the pendulum arm. When the pendulum bal¬ 
anced, the equation of moment of force was 

Gj OAI sin a l = F/.IOM L | cosa L (1) 

Then 

FjOM L l= GjOAltancx L (2) 

As tree model fixed, the pendulum system was 
forced by 

(1) Gravity of pendulum arm G L \ 

(2) Drag of pendulum arm Fl. It’s center was at the 
point Ml of the arm. The influence of the attached 
tree model on the flow forcing on pendulum arm was 
so small and the difference between angle a L and 
a T were also so small that follow equations could be 
truth approximately 


Fl= F l 

(3) 

Iom l I = IoMl’I 

(4) 

and 
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F l I OM l I = F l ’ I OM l ' I = G l I OA I tan a L (5) 


o 



(3) Gravity of tree model Gy. Its center was at middle 
point M 0 of the model. 

o 



Fig. 1. Force diagram of the pendulum system (Swing angles were exaggerated ) 


(4) Drag of tree model F T . Its center M T was related 
to the vertical distribution of permeability of the tree. If 
the permeability was uniformly distributed, the height 
of point My could be estimated by follow (in neutral 
stratification) 



In —) 2 adz = 
z o 


ecu. 

2k 1 


H 2 H . 2 

(In —Y -2ln—+2- 

z 0 z o H 


( 6 ) 


where u *, z 0 and * were friction velocity of ap¬ 
proaching flow, roughness length of the floor and 
Karman constant, so IOMtI=2IOAI+(H-Ht). 

The balance equation of moment of force was: 

Gi IOAI sin a r*-Gy(2 IOA l+0.5H)sin a T = 

Fl’IOMl’ I cos a t+ Ft lOMy I cos a y (7) 


Error analyses of the experiment 

Errors could be caused probably from follow: 

(1) Displacement reading: the longer as the dis¬ 
placement was, the bigger the reading error was. It 
was generally about 1 mm. The average displace¬ 
ment was about 4 cm and the relative average read¬ 
ing error was about 2.5%. 

(2) Suspending of the tree model: the gap between 
trunk base and the tunnel floor was about 1 mm. This 
corresponded to omit 1 mm high trunk base. Wind 
velocity here was very low and the omission impacted 
little of drag (only 10 4 of the drag). So the suspending 
errors could be neglected. 

(3) Dip and elevation of the tree model: the maxi¬ 
mum elevation of the model in the experiment was 
0.15 cm, corresponding to 2% of the model height. 
This caused the errors of drag about 2%. The maxi¬ 
mum dip angle of the model was 2.6°. This may be 
neglected comparing with the dip of actual tree. 

Results and discussion 


Considering (5), then 

Fj- (Gi IOAI sin a y+Gy(2 IOA I+0.5H) sin cty — 

Gl IOAI tan a L cos a T )/( lOMyl cos a y) (8) 

All the terms in the right hand were measurable. So 
the drag of tree model Fy was easy to calculate. 


Drag of isolated tree model 

The measuring result showed in Fig.2. It shown that 
drag of isolated tree model was related to its perme¬ 
ability (or porosity) and approaching wind velocity. 
Less permeability and stranger wind velocity gener¬ 
ated bigger drag. But the drag of cylinder model was 
less than that of expected. 









GUAN De-xin et al. 


159 



Fig. 2. Variation of drag of isolated tree Frwith 
permeability a under different approaching wind 
velocity at tree height 

Drag coefficient of isolated tree model 

Drag coefficient of isolated tree mode was defined 
as 

C t = 2F t I (DHp u e 2 ) (9) 

where Ft and C r were drag and drag coefficient of 

isolated tree respectively. D and H were the crown 
width and height of the tree, p was air density. u e 
was geometrical mean wind velocity of the flow 


1 H r u. z , _ u» . H 

u e — ——■ I In dz In— (10) 

H *' K V ot 


K 


ez n 


Fig.3 showed the variation of drag coefficient of 



Fig. 3. Relationship between drag coefficient Crand 
permeability a of isolated tree (dots were measure¬ 
ments and line was the empirical model) 

isolated tree versus its permeability. Except for cyl¬ 
inder model, the drag increased with decrease of the 
permeability. Cylinder model had a slight small drag 
coefficient compared to dense tree model. This were 
caused by the attribute of cylinder model. It had 
smoother surface than that of tree model. And the 
friction with the flow was less. The relationship be¬ 


tween drag coefficient of isolated tree and its per¬ 
meability was empirically described by equation 

C T = -1.251 a T 2 + 0.489 a T +0.803 (11) 

with correlation R^O.981. 

Conclusions 

As the single tree was the basic element of inter¬ 
action between atmosphere and some representative 
surface such as Savannah, agroforest and so on, 
drag of isolated tree was an important index for the 
aerodynamic research of these surfaces. 

Because no instruments were available for meas¬ 
uring so small drag of the tree models in wind tunnel. 
The pendulum method was designed, and drag of the 
isolated tree model could be calculated by the for¬ 
mula evaluated from balance equation of moment of 
force. The experiment was easy to conduct and the 
errors were small. 

The results showed that besides the crown scale, 
drag of isolated tree increased with decrease of its 
permeability or porosity. The drag increased as ap¬ 
proaching wind became stronger. Drag coefficient of 
isolated tree decreased as its permeability increased. 
This relation empirically described by quadric curve. 

This paper gave the results obtained from wind 
tunnel experiments. The tree models were relatively 
ideally constructed and in some extent, were different 
from actual trees in geometry and physical attributes. 
So the drag measurements of actual trees in field are 
expected. Nevertheless, the results in this paper 
gave the basic law about drag of isolated tree and will 
helpful for further researches. 
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